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The correlated two-particle problem is solved analytically in the presence of a hnite cavity. The 
method is demonstrated here in terms of exactly solvable models for both the cavity as well as 
the two-particle correlation where the two-particle potential is chosen in separable form. The two- 
particle phase shift is calculated and compared to the single-particle one. The two-particle bound 
state behavior is discussed and the influence of the cavity on the binding properties is calculated. 
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INTRODUCTION 

Finite quantum structures are intensively studied con- 
cerning their transport properties, trapping behavior as 
well as considering them as artificial atoms 0, 0| which 
become especially interesting in reduced dimensions. The 
main tool to investigate the properties of such quantum 
cavities is the scattering of particles connected with ex- 
citations || or the current measurement [1.0.0] through 
these structures. We want to regard as quantum cavities 
here also quantum dots or any finite structure. The ef- 
fort of theoretical investigations has reached a fairly high 
level adequate to the obtained experimental quality of 
data. Nevertheless it is noticeable that almost all inves- 
tigations are restricted to one-particle transport proper- 
ties like current, conductance etc. On the other side the 
experimental abilities are on such level that also higher- 
order correlations and more than one-particle properties 
can be envisaged. The aim of the present paper is to 
investigate two-particle correlations together with one- 
particle correlations on the same footing. Since these 
investigations suggest to look at a new type of experi- 




ments on quantum cavities we want to present an ex- 
actly solvable model for exploratory reasons. In this way 
the method is demonstrated which can then be used for 
any realistic nanostructure, with an amorphous or crys- 
talline system as background provided the one-particle 
wave functions are available from e.g. standard ab-initio 
methods. 

As illustrated in figure ^ we consider the gedanken ex- 
periment to scatter two particles with each other in the 
presence of a quantum cavity. In usual approaches it 
would be sufficient to know the scattering properties of 
single particles (though correlated) on the cavity in order 
to describe the transport properties. We want to demon- 
strate that the two-particle correlations in addition to 
the scattering on the quantum cavity lead to interesting 
new effects of coupled channels which might be worth an 
experimental investigation. 

As sketched by the figure, we face a two-center scatter- 
ing which is conveniently described by coupled-channel 
scattering theory Q. A similar problem has been con- 
sidered by Beregi et al. who solved the pro blem with 
the help of the Faddeev technique 0, 0, . We will 
adopt here a simpler approach which is possible due to 
the Gell-Mann and Goldberger formulation Due to 
the different occurring interactions, we define two chan- 
nels, the one-particle scattering channel V\ = V a + Vb and 
the corresponding two-particle one, V2 = V a b- The total 
T- matrix can be found by the Gell-Mann and Goldberger 
two-potential formula summarized in appendix „ 



T = T 1 + (1 + V 1 Q 1 )T ab {l + GM) 
Here the one-particle T-matrix is given by 

T 1 = v 1 + v 1 g 1 v 1 



(i) 



(2) 



FIG. 1: Schematic picture of correlated two-particle scatter- 
ing on a quantum cavity. 



and the correlated two-particle T-matrix T ao is a solution 
of the modified Lippmann-Schwinger equation 



Tab = V 2 + V 2 GiT ab 



(3) 
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where the full solution of the one-particle channel enters 
in Q\ . For that one assumes that the one-particle problem 
is solved and the propagator in the one-particle channel 
is known Qi = (uj — Hi)^ 1 . This scheme is advisable if 
one channel can be solved exactly and the other channel 
is then described by a Lippmann-Schwinger equation © 
with the help of the solved channel. This method has 
been compared to other coupled-channel methods [l2T ] 
and successfully applied to defects in potentials 0] and 
in developing a hybrid representation of the nuclear po- 
tential 14] . We will use this scheme here and will employ 
an exactly solvable model for both channels. 



MODEL 
Model cavity 

For the cavity we choose a model of an opaque wall 
in three dimensions which can be found in text books 
[lfij . This model is able to simulate properties of opaque 
confining traps as well as parametric resonances. The 
potential in the single-particle channel takes the form 



V a (r) = ~8(r-R) 
2m R 



(4) 



with the reduced mass m of particles a and b scaled out. 
Thus we have an opaque delta potential at the radial 
distance r — R from the center of cavity. The wave 
functions out- and inside the cavity decouple completely 
in the limit of hard spheres, ft — > oo, where the interior of 
the cavity is characterized by discrete energy levels. For 
the opaque wall, ft < oo, outside and inside are coupled 
and resonances occur at the virtual levels. 

The wave function of this radial problem has the form 
$„(r) = (r\n) = Xn( r )/ r - While the subscript n denotes 
general quantum numbers the specific wave function here 
shall be characterized by continuous momenta n — k such 
that 



Xk(r) = c 



A kR sin kr 



;r < R 
;r > R 



sin (kr + . 

The one-particle phase shift <p is determined by 

cot (x + (j>) — cot x = ft/ x, 
where x = kR. The amplitude can be written as 



2 _ sin (x + , 



1 + ta,n z x 



■ 2 

sin x 



tan 2 x + (1 + ^t&nx) 2 
and shows resonances at the virtual levels 



tan 2x,, 



(5) 



(6) 



(7) 



(8) 



with even v which approach the levels x„ = k v R = vix 
in case of an impenetrable wall, ft — > oo. Assuming only 




FIG. 2: The single-particle cross section on the cavity JjJJ with 
zeros at kR — nn (solid line) and the amplitude Q (dashed 
line) versus momenta for = 20. The filled circles marks the 
maxima of the amplitude (Sf, the open circles the maxima of 
the cross section according to IjlOjl and the open squares the 
minima and maxima of the cross section according to 1111 . 



s-wave scattering the corresponding single-particle cross 
section reads 



a(k) = 



4tt 



k 2 1 + cot 2 i 



(9) 



which has zeros at x n = nn. From © one finds that the 
cross section has a limited number of maxima at 



2— = sin 2x u 
ft p 



(10) 



depending on the potential strength, ft — 7r/2 > 27ij/i|. 
These maxima are accompanied by minima and an infi- 
nite number of further maxima at 



— 2— = tan 2x u . 
ft M 

These features are illustrated in figure |21 



(11) 



Two-particle problem with separable potential 

We study in this paper the combined problem of single- 
and two-particle scattering and will use a separable rep- 
resentation of the two-particle interaction which provides 
tremendous technical simplifications. One can show that 
the effect of any finite-range potential on the scattering 
wave function \I/ can be represented in separable form 

mm, 



(12) 
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where \p) is the undisturbed plane wave and gi are form 
factors. This is known as the expansion of the wave func- 
tion inside the potential range and outlined in appendix B 
Separable potentials are powerful tools for constructing 
unknown nuclear matter potentials from available exper- 
imental scattering phase shifts by inverse scattering the- 
ory [l^. The application of such separable potentials 
is versatile ranging from plasma physics [lflj ] to medium 
modification of transport properties in nuclear matter 
Efllilljfiij l up to high field problems in superconductiv- 

According to the Gell-Mann and Goldberger formula 
(|T]l we assume the one-particle wave function \ni) to be 
known 



Hi\ni) = E ni \m). 



(13) 



The two-particle T-matrix equation J3J is now repre- 
sented in momentum states. We note the momentum 
difference of the two particles by small letters, e.g. p, and 
the center-of-mass momentum by capital letters, e.g. P. 
Introducing a complete set \n1n2) according to l|13|l we 
obtain 

( Pl P\T ab \ P2 P') = V P1P2 5 PPI 

ST t/ (pPln^jn^lp'P) , ,p |/r . 



n\n 2 pp' P 



UI — E n 



11] 



(14) 



By convention we use — I (if^p" an< ^ set H= 1 within 

this paper. Equation (|14fl is the exact formulation of the 
quantum two-particle problem in channel 2. It can be 
solved with the help of such a separable potential [24| 



Vpq = ^9p9q 



(15) 



with form factors g p as given in appendix B While gen- 
erally any finite-range potential can be represented in a 
separable form of finite rank, we restrict ourselves here 
to rank-one potentials for exploratory reasons. The tech- 
nique remains the same for high-rank potentials. 



TWO-PARTICLE SCATTERING ON FINITE 
CAVITY 

Free two-particle problem 

The case without cavity or unrestricted two-particle 
scattering is solved by plane waves with the center-of- 
mass momentum K and the difference momentum k as 
good quantum numbers, \n1n2) = \kK). Performing the 
k, K integrations, the Lippmann-Schwinger equation l|14|) 
reads with (pxP^a^P') = Spp'tp 1 p 2 (P) 



tpip 2 (P) — ^9pi9p2 + ^ ^ 9pi 



t. 



9q- 



<U>1 



OJ - 



P 2 V 1 
x — + 1>V 

4m m ' 



(16) 



which is easily solved as 



^PlP2 (P) 



1 - XJ(P,lu) 
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p 2 <? 2 1 • ■ 

OJ — T 1 IT) 

4m m ' 



(17) 



This result is given explicitly for different form factors in 
appendix 



Two-particle scattering and cavity 

Now we are going to consider the disturbed two- 
particle scattering problem, which should reduce to 1|17|) 
in the limit of unrestricted two-particle scattering. The 
equation for the T-matrix (|14f> can be simplified by in- 
troducing the modified form factors 



P 

Splitting off the trivial factors 

( P iP\T ab \p 2 P') = \g pl g P2 t(PP'), 



(18) 



(19) 



from (|14fl the equation for the center-of-mass dependence 
follows 



t(PP') = 5 PP , + A £ 



G„ 1 „ 2 (P)EG; in2 ( J P)i(PP') 



w - E ni7l2 + irj 



(20) 



The completeness of \n1n2) implies the relation 

£ G nin2 {P)G* nin2 {P') = 5 PP ,J29 2 q - (21) 

n±n 2 q 

We introduce the auxiliary function 

K 1 n 2 (P')=Y, G n 1 n 2 (P)t(PP') (22) 
P 

from which the T-matrix can be obtained with the help 
of (EU as 

t(PP') = -L_J2 G nrni(P)K in2 (P'). (23) 

Multiplying (|2T]|l with G* nim2 (P) and integrating over P 
we obtain an equation for the auxiliary function <|22|) 

^ l m 1 m 2 (P ) = G mi „ l2 (P ) + Hm 1 m 2 n 1 n 2 h nin2 (P ) 

n±n 2 

(24) 
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with 



(P) 

P 



iii i in ~2 n i n-2 



(25) 



For cavities possessing only discrete states in which 
case n, and are natural numbers, (|24|l is a linear equa- 
tion system which can be easily solved. With the help of 
(|23[l and (|19fl the T-matrix is consequently found exactly. 

The situation becomes more complicated if continuum 
states are present in the cavity. Then the quantum num- 
bers {n\n2\ becomes continuous variables like in the case 
of the free two-particle problem and l|24(l turns into an 
integral equation. Wc cannot present at the moment the 
analytical exact solution of the combined problem con- 
sisting of bound and scattering states. The numerical 
scheme consists in discretizing this continous integral and 
solving the linear equation system (|24|l . This complete 
numerical solution will be devoted to a forthcoming work. 
To get a first physical insight, we will instead restrict 
ourselves in this paper only to cavities with no bound 
states, i.e. only scattering states, and moreover to scat- 
tering events where the center-of-mass momenta of the 
two particles are not remarkably changed. 



Approximated solution 

Due to the spatial structure of the one-particle prob- 
lem, in general we will not have translational symme- 
try. This means that the two-particle problem need not 
to have the same center-of-mass momenta of incoming 
and outgoing particles. Only in special cases like the 
zero center-of-mass momentum and no change of the to- 
tal momentum due to release or capturing of momentum 
by the cavity, we can assume approximately a diagonal 
form t(PP') w Sppd(P). 

We understand this approximation as considering the 
effect of the cavity on the free two-particle problem l|16(l 
and (fTTj) in first order. Such a diagonal solution is ob- 
tained, if for the momentum integral over P on the right 
hand side of (|20(l we apply the mean value theorem of in- 
tegrals and approximate the mean value P by the center- 
of-mass momentum P 



"£G* nin2 (P)t(PP') 
p 



t(PP')5>* lTi2 



(P) 



t{PP')^Gl 
p 



,(P). (26) 



Then H2U|) is easily solved and the diagonal T-matrix is 

(27) 



with 



j(p» = J2 



J2 G ni n 2 (P)G* nin2 (P') 



E ni n 2 



IT) 



(28) 



Alternatively we can obtain this result by assuming 
ad-hoc a diagonal T-matrix t(PP') w 8ppd(P) which 
according to f2T)|l fulfills the integral equation 



t{P')5 PP , = 



)p P , 



G ni „ 2 (P)EG; in2 (p)i(p')^P' 



W - E nrn2 + U] 



(29) 



Integrating this equation over P we obtain the same ap- 
proximate solution (12711 



t(P') 



1 - \J(P',u)' 



(30) 



It is interesting to note here that (|30|) is also a solution 
of a nontrivial integral equation. To see this one can 
integrate H29[) over P' instead over P as done above. One 
then obtains the integral equation 



t(P) = 1 + A£ 



G ni n 2 (p)EG; lin2 (P)t(P) 
p 

w - E nifl2 + if] 



(31) 



Therefore we can conclude that ptljl is also an approx- 
imate diagonal solution of the integral equation 
which has not been known previously to the authors' 
knowledge. 

The approximate solution Ij27(l together with (|28|l and 
<P)l is the main result which we will employ within this 
paper. It represents a generalization of the standard sep- 
arable T-matrix result, e.g. of Ref. |24|, in two respects: 
(i) It describes the two-particle scattering on a finite cav- 
ity or on an arbitrary structure represented by the one- 
particle wave functions and (ii) it gives the exact solution 
of the two-particle quantum problem in the case of con- 
servation of total momentum before and after scattering. 

Before continuing we want to show first how the limit 
of unrestricted two-particle scattering is reached. Unre- 
stricted scattering means the absence of the cavity so that 
the single-particle wave functions become plane waves, 
\n1n2} = \kK) as discussed in section D Consequently, 
(|18fl takes the form 



G^(P) = g k Sp K 



and J5HJ is 



J inf (P W ) = £, 



Ei 



+q 



- E, 



(32) 



(33) 



which is the standard expression for unrestricted two- 
particle scattering 117|) and separable interaction. In nu- 
clear physics this case is labelled infinite matter. 
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ANALYTICAL RESULTS FOR TWO-PARTICLE 
SCATTERING ON AN OPAQUE WALL 

We are now going to evaluate the result l|77|l with 
and (|18fl . Restricting the calculation to a rank-one sep- 
arable potential characterized by only one form factor, 
e.g. g p = exp (— p 2 /4/3 2 )//3 2 , and the potential strength 
A, one can describe two properties of the two-particle cor- 
relations, the scattering length ao and the range of the 
potential ro as illustrated in the appendix B For clarity 
we will first concentrate only on the scattering length and 
will separate the effects of potential range. Therefore it 
is advisable to start with the contact potential which is 
represented by the limit of zero range (3 — > oo. In a sec- 
ond step we will discuss the effects of finite range of the 
potential. 



Contact potential for two-particle correlations 

The separable potential allows us to adopt the contact 
potential by a special limit using a parameter-dependent 
coupling constant A(/3) such that the scattering length 
is kept fixed. In this way one ensures intrinsically the 
necessary renormalizations which otherwise have to be 
performed by subtracting diverging terms to fix the scat- 
tering length to the physical value |25| . 

Mathematically, the limit (3 — > oo with fixed scattering 
length is realized if we rewrite the coupling constant A in 
terms of the scattering length. This is obtained from the 
general relation between phase shift ip and the T-matrix 



tan <p 



ImT 



RcT 



inf 



A4 nf 



(34) 



where we will abbreviate the real and imaginary parts of 
the unrestricted two-particle scattering expression i|33|) 
as J^ lf = ReJ inf (P, and J{ ni = ImJ inf (F, £\ respec- 
tively. One knows that for contact potentials 26] 



tan ip" 



pa 



(35) 



for all momenta p and we express the parameter- 
dependent coupling constant A(/3) from (|3*4l as 



A 



pa 



J{ nt 



Inserting this expression into 127J1 we obtain 

^9pi9 P 2 _ pao9 Pl g P2 



(36) 



T ab (P,iu) = 



1 - \J(P,U>) 



Tinf 
J I 



pa a (J^- J(P,w))' 
(37) 



Now we are ready to perform the limit of zero range 
(3 — > oo. From l|28|> and (|33|) one sees that the form 
factors g p = l/(/3 2 + p 2 ) appear with the same order 



in the numerator as well as in the denominator which 
cancel in the limit of zero range. Consequently, (|37|l for 
the T-matrix l|27|l approaches 



T ab (P,uj) 



pa 



J? inl - pa ( Jg inf - J°(Poj)) 



(38) 



where the additional superscript denotes that both ex- 
pressions H28JI and H33[) have to be taken with g p — > 1. 
The difference J^ nl — J° in the denominator Q38J1 ensures 
the convergence of the expression since Re J and Jp lnf di- 
verge. 

It is instructive to check the T-matrix for contact po- 
tentials without cavity. From (|33J) one has 



p 2 



i 



(/>+ 4^ - 2^ 



dkk 2 

UJ - — 



(39) 



which leads to the T-matrix l|38|l for contact potentials 
in unrestricted two-particle scattering 



1 



-— + ip 

a r 



4tt 



VmEb ; uj = -Eh 



(40) 



Though (|39|l diverges, the solution l|38|l remains finite. 
With the help of (|34|) one can convince oneself that l|35|) 
is fulfilled for positive energies uj — p 2 /m as it should 
be and that in the case that the scattering length 
is negative indicating attraction, the bound state E b = 
l/ma 2 , represents the pole of the T-matrix at uj = —Eh- 
Now we are going to calculate J° according to in- 
cluding the finite cavity represented by the single-particle 
wave function \n) in JSJ). The explicit form for l|18|l reads 

G°n in2 (P) = $>P|nm 2 >= f dre- irP * ni (r)*„ 2 (r) 
p J 

(41) 

where the spatial representation of the one-particle wave 
function is ^ ni (r) = (r|ni). This means that 

J2 G °nn 2 ( P ) = ^(0)^(0) 



G° 1K2 (0) = V25 ni1 



(42) 



The last equality follows from the completeness of the set 
of single-particle wave functions as well as the normaliza- 
tion 



J dr\V ni \ 2 = V2 



(43) 
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ensuring that the normalization of the wave function 
\n\n2) is 2 for the two particles. 

Now we can give a simple expression for (|28f) in the 
case of center-of-mass scattering P = which we want 
to discuss in the following. We obtain with the help of 

62J 



J°(o,-) = £ 



l^(0)| 2 V2 



2E„ 



IT] 



(44) 



Together with l(39() we have obtained the two-particle 
T-matrix for contact interaction ((38(1 provided the one- 
particle wave function of the cavity ^ ni (r) = (r\ni) = 
Xni( r )/ r an( i the single-particle energy E ni according to 
((13(1 are known. For this purpose we use JHJ) of the opaque 
wall and the necessary part of the wave function in 144(1 
reads 



* ni (0) = ckA 



kit 



(45) 



with the amplitude given by 10. We determine the nor- 
malization c in such a way that the particle should be in 
a sphere with radius L which is set to infinity in the end. 
For such large boundaries we can employ radial cyclic 
boundary conditions, xi r + L) = x( r )i such that we have 



dk. 



(46) 



Calculating the norm 1(43|) with JSJ it follows that 
Lc 2 72 



lim 

L — >oo 2,7T 



At: 2 ' 



(47) 



Therefore, the arbitrary large fictitious boundary length 
L is dropping out in 144(1 and we obtain 



oc 

Re j°( ,„) = ±J 



k 2 A 2 
dk k AkR 



1.2 

LU - — 



ImJ>,^ = ^) = -^A 2 



An 



pit- 



(48) 



This allows us to determine the T-matrix ((38(1 with the 
help of as 



z= L ^(oA) = (-i-^ ( P fl) 

Anao m V R 

with 



ipa A 2 



pit 



dx- 



(49) 



(50) 



Please note that the —1 term comes from the subtrac- 
tion of ((39|) in ((38(1 which renders the expression finite. 
The function F can be given analytically as presented in 




FIG. 3: The one-particle phase shift (top) of an opaque wall 
@ versus wave vector and the free two-particle phase shift 
135H as well as the combined ones 15111 . The wall parameter 
is O = 10 which corresponds to Jh(0) = 0.91 and the two- 
particle scattering length ao = ±0.LR. In the bottom figure 
the corresponding cross sections 152H are plotted. 



appendix , 
shift 



From ((49(1 we obtain the two-particle phase 



tan<^(p) 



pR A 2 



pit 



F n (pR) 



(51) 



In figure we plot the phase shifts for a repulsive 
opaque wall with f2 = 10 and a weak two-particle scatter- 
ing length of ao = ±0.1i?. The one-particle phase shift cf> 
is negative indicating repulsion and has the typical oscil- 
latory structure [l5T | due to the resonances at the virtual 
levels. The free two-particle phase shift ip mi of 1(35(1 ap- 
proaches 7r/2 for large k and starts from zero for positive 
ao and from tt for negative scattering length indica ting 
a bound state according to the Levinson theorem [28j. 
The two-particle phase shift 1(51(1 for the case including 
the opaque wall follows the free two-particle scattering 
case for such small scattering lengths though the oscilla- 
tions of the one-particle problem are visible. Interestingly 
these oscillations are shifted slightly when positive and 
negative scattering lengths are compared. 

With the help of ((51(1 the corresponding cross section 
takes the form 



An 



P 



j sin <p 



AtxR' 



- 1 ,' 



x*A* + (F n (x) + %? 



(52) 



which will be determined by the structure of the one- 
particle amplitude A and a combination of A and F de- 
pendent on ao in the denominator of 1(52(1 . The latter 
leads to additional structure in the cross section and can 
be understood as an interference between the one-particle 
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and two-particle properties. This can barely be recog- 
nized in figure |31 where a second maximum appears be- 
sides the oscillations given by the one-particle cross sec- 
tion. The appearance of this second maximum is much 
more pronounced when we go to larger scattering lengths 
corresponding to stronger two-particle correlation as in 
figure 0] Therefore, it is a correlation effect in the two- 
particle channel. 

Figure 0] also shows that the two-particle problem in 
the presence of a cavity has no bound state anymore 
though the free two-particle problem has one. The re- 
pulsive opaque wall has destroyed the bound state forma- 
tion. In order to understand this behavior we consider 
the T-matrix not at positive energies to — k 2 /m as done 
so far but at negative energies where the pole yields the 
bound state energy uj = —Ef,. The T-matrix l|38|l reads 
in this case 



-Airao/m 



-(53) 



1 + a yfmE h + f FaiiRy/mE] 

with F from i|50[l . Consequently, the bound state energies 
are given by the pole 



— = F n (ib) + b 
a 



(54) 



with b — Ry/mEi,. It is easy to prove that Fq (*i>) + b > 
Fq (0) , the explicit analytical expression can be found in 
appendix B This means, as soon as the cavity parameter 
17, R and the two-particle scattering length a fulfill the 
inequality 



R 

«o 



> F n (0) 



(55) 



we will have two-particle bound states. In figure this 
is the case for negative scattering lengths since R/clq = 



— 10 and -Fio(O) = 10/11 while no bound states occur 
for positive scattering length R/clq = 10. Oppositely, in 
figure 01 for both positive and negative scattering lengths 
we have R/ao = ±0.1 < .Fio(O) = 10/11 and no bound 
states appear anymore. In this figure the free two-particle 
problem has a bound state for negative scattering length 
while the repulsive wall allows no bound state for these 
parameters. According to 1)55(1 the limiting scattering 
length is ao = — 1.1R for £1 = 10. In figure El we plot the 
cases just slightly above and below that critical value. 
One sees that the bound state present in the free case is 
destroyed by the cavity if the scattering length is more 
negative than the lower critical value. 

It is now interesting to consider the attractive wall fl < 
0. Then we might expect that bound states are favored 
in any case which is surprisingly not true. 

We plot the function Fn(0) in figure HJ In the range 
ill < < 0, -Fh(O) becomes negative. Therefore any 
negative ao leads to a bound state according to l(55)) . 
Both two-particle scattering as well as the opaque wall 
are attractive and also their combination acts attrac- 
tively. Positive do will lead to a bound state if ao > 
— R/Fq(0) > 0. Though the wall is attractive to both 
particles, bound states occur only if the two-particle cor- 
relations are larger than a critical value. 

This is illustrated in figure In this figure we have 
chosen free two-particle scattering parameters which do 
not lead to bound states. The attractive opaque wall 
creates a bound state for ao = 0.5i? while no bound state 
occurs for a = 0.25R. 

For the range f2 < f2i the function .Fh(O) is positive. 
This means that the same discussion holds as for pos- 
itive fl. According to l(55|l ao must be negative and 
ao > —R/Fq(0). This means that bound states occur 
only if the two-particle scattering length is larger than a 
critical value. The appearance of such critical values for 
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FIG. 6: The function .Fh(O) of 150I I which allows bound states 
if —R/ao > Fn{0) according to l|55fl . 



the two-particle scattering length is an unexpected effect 
of interference. Please note also that if 1 < Fq(0) < 1.21 
there are two values of £1 which lead to the same Fn(0). 
This means we can have two different cavities which 
yield the same bound state. 



8 single particle £2=-1 .2 
a =0.25R, £1=0 
a =0.25R, £2=-1.2 
a =0.5R, £2=0 
- a„=0.5R, £2=-1.2 




FIG. 7: Same as figure El but with Q — —1.2 which corre- 
sponds to Fq(0) — —3.33 and an — 0.25R and an = 0.5R. 



Finite range effects 

After discussing the effect of the scattering length we 
want to return to the more realistic situation of finite- 
range two-particle interaction. Therefore we choose a 
separable form of the potential and evaluate (|28|) explic- 



itly. For exploratory reasons we choose the exponential 
form factor of appendix D and the center-of-mass system 
P = 0. 

Introducing the dimensionless integrals 



Gnin 2 (0) : 



— 3 ffa^kt (q) w k2 (— q) = —p-G {xi,x 2 ) 



(27T) 



dp 32pRc* - 
— ^G niri2 (p) = 7 ^G(xi,x 2 ) 



with x n = Rk n , we obtain with the help of H46() and 147|l 

^ T(n *\ 2 6 £/? 2 i? 2 L a G (x 1 ,x 2 )G(x 1 ,x 2 ) 
AJ ( ' W )= T 5/ 2 I<bidx2 _ „ 2 _ „ 2 , jz -( 57 ) 



(56) 



LO — XT — X% 



i q 



Here we used lo* = 2mR 2 Lo and the dimensionless cou- 
pling constant £ which can be written in terms of the 
scattering length ao according to iffify 



mX 



1 



(58) 



Consequently, the potential is attractive for scattering 
length ao > or ao < — ^J-k/2/ (3 and repulsive otherwise. 
The potential range according to (|86|l follows from l|93|) 

as 



_1 I 2 1 



(59) 



which shows that for negative scattering length we should 
have ao < —^/tt/2/2/3 to render the potential range pos- 
itive. 

It is convenient to use the spectral representation fol- 
lowing from Ij57(l 



ReAJ(0,w* 



1 f , ImAJ(0,w) 
— / duo 

7T J LO* — LO 





lm\J(0,Lo*) = - 



2 5 £0 2 R 2 



r 3/2 



x / dxi 

o 



Gq(xi, y/w*-x\)G{X\, y/b)*-xl) 



^Jlu* — x\ 



(60) 



The explicit analytical integrals Go and G can be found 
in appendix B 

The results are presented in figures l%l 101 We compare 
the case with a cavity to the free two-particle scatter- 
ing. The limit of contact potential (3 — > oo of the last 
chapter is plotted for comparison too. In figure [S] we 
have chosen a scattering length which does not allow the 
formation of bound states. The two-particle phase shift 
follows the general trend of the phase shift without cav- 
ity but shows characteristic oscillations which are visible 
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FIG. 8: The two-particle phase shift (top) and cross section 
(bottom) in the presence of the cavity (thick lines) and with- 
out cavity (thin lines). The contact-potential case (3 — > 00 
(dotted lines) is compared with the finite-range case (solid 
and dashed lines). The scattering length is ao = R. 



3.0 




kR 



FIG. 9: Same as figure|Hl but for ao = -4R. 

in the cross section too. In the figure we chose such a 
scattering length that the free two-particle scattering has 
a bound state. While for the contact potential the repul- 
sive opaque wall already destroys the bound state, the 
bound state survives in the finite-range potential. This 
means that the finite range enhances the correlations. 

Decreasing the potential strength of the attractive two- 
particle scattering by a more negative scattering length 
the repulsive opaque wall will destroy also the bound 
state for the finite-range potential. In figure ITU1 we plot 
the limiting case where the bound state is just still real- 
ized. 

From the figures it is obvious that the finite-range po- 



FIG. 10: Same as figureEl but for a = -WR. 

tentials are more realistic than the contact potentials. 
For higher energies the phase shifts for the contact po- 
tential approach an artificial value of 7t/2 while the phase 
shifts for the finite-range potential approach zero as they 
should. For finite-range potentials also the cross section 
falls off faster than for contact potentials. 

SUMMARY 

We have investigated the coupled-channel problem of 
two particles interacting with themselves and with a cav- 
ity. We have chosen exactly solvable models for both the 
cavity as well as the two-particle potential in order to 
demonstrate the method. With the help of the single- 
particle wave function of the cavity, the two-particle 
problem is solved. For cavities with bound states the 
T-matrix is determined exactly. In the case of additional 
scattering states we restrict ourselves only to processes 
where the change of the center-of-mass momenta are neg- 
ligible and present an approximate solution. 

The influence of the cavity on the two-particle prop- 
erties like scattering phase shift and bound states is in- 
vestigated. We observe oscillations in the single-particle 
and two-particle phase shifts. While the oscillations in 
the single-particle phase shift are the standard resonances 
at the virtual levels of the cavity, the two-particle phase 
shift shows besides these oscillations also interferences 
which are visible in a splitting of the maxima in the cross 
section. 

For repulsive opaque walls the scattering length must 
be negative but larger than a critical value in order that 
bound states remain. For attractive walls there is a spe- 
cific range of negative wall parameters for which any neg- 
ative scattering length will lead to bound states since 
both the cavity wall and the scattering length act attrac- 
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tively. In addition also positive scattering lengths will 
lead to bound states in this case provided the scattering 
length is larger than a critical value. For more attractive 
walls we could show that bound states are only possible 
for negative scattering lengths which are, however, lim- 
ited by a negative critical value. The occurrence of such 
lower critical values is an unexpected effect of interfer- 
ence. 

The model discussed here has the merit to be analyt- 
ically solvable in all parts and serves as a toy model for 
more realistic situations. In particular, the finite cavity 
can be easily described by a more realistic wave function 
which enters the two-particle problem in such a way as 
worked out here. We intended to provide in this way a 
clear method which is also easily applicable to the case 
of two-particle scattering on any finite or more realistic 
structure, like crystals or amorphous structures. In this 
case one has to provide the single-particle wave functions 
e.g. by an ab-initio method and can use the method here 
to calculate the two-particle correlations explicitly. 

The discussion here has been limited to the two- 
particle correlations for exploratory reasons. In a forth- 
coming work we will show how this method is applicable 
for the situation of many-particle correlations. 

This work has been supported by DA AD. The discus- 
sions with Ingrid Rotter which have been of great value 
are gratefully acknowledged. 



V = and can be written as a sum of auxiliary T- 

i 

matrices 

TH = ]T Vi + X>i0o(a;)TH = E^H' ( 65 ) 

i i i 

These auxiliary T-matrices are defined by 

T{ = Vi + ViQoT (66) 

and can be rewritten as 

T( = Vi + ViQiT-ViQiVigoT 

= Vi + ViQiVi^ViGiVjil + QoT) 

= Ti + J^ViGiTj (67) 

where we used (|63[) to write the first line and the defi- 
nition H66|) to derive the last equality. Using further the 
identity ViGi = %Go the final form reads 

7>) = 7^)+EWSoWH (68) 

which shows how the single particle channel T-matrices 
% are coupled to yield the auxiliary T-matrices which 
add up to the total T-matrix (|65|l . 



FORMULAS FROM COUPLED-CHANNEL 
SCATTERING THEORY 

Many coupled channels 

We consider a system consisting of i channels. Using 
the free propagator 



1 



uj - Ho + it) 



(61) 



and the propagator including the interaction of the zth 
channel 



GM 



1 



oj — Hi + ii] 



(62) 



the Lippmann-Schwinger equation for the ith-channel 
propagator reads 

Qi(u) = Q (u) + Gi(u)ViGo(u) (63) 

corresponding to the ith-channel T-matrix 

TM =Vi + ViQi(oS)Vi. (64) 

The total T-matrix is described by the Lippmann- 
Schwinger equation with respect to the total potential 



Derivation of the Gell-Mann and Goldberger 
formula 



In the case of only two channels the coupled T-matrices 
read 



and 



T = T/ + T 2 ' 
T/ = T + TxGqT 2 
% = T 2 + T 2 g T{ 



T = Ui + ViSoTi 
T 2 = V 2 + V 2 GoT 2 . 



(69) 



(70) 



Introducing T/ into T 2 from l|69(l leads to 

T 2 ' - T 2 (l + goT 1 ) + T 2 g T 1 g T 2 ' (71) 

with the help of which we define 

T ab = T 2 \l + Go^)- 1 = T 2 + T 2 GoT 1 G T ab . (72) 

The total T-matrix from l|69() can then be written as 

T = T 1 + (l + T 1 g )T 2 ' 

= T 1 + (l + T 1 g )T ab (l + ^ T 1 ). (73) 
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It remains to show that T a b is the T-matrix obeying 
a Lippman-Schwinger equation with the propagator Gi- 
For this purpose we use ijfiSl) 

— Go = —Q\ + GoViGi 
{l-V 2 Gi)- 1 V 2 = (1 - TzGoV^i)- 1 ^ (74) 

where the second line follows from the first by trivial alge- 
bra. Now the right hand side of l|74l) is just the definition 
of T ab in |Z3) using once more TGo = ViGi- Therefore 
we see that from l|74|) the Lippmann-Schwinger equation 
© for Tab follows, 



T ab = V 2 + V 2 GlT ab 



(75) 



which completes the proof of the Gell-Mann and Gold- 
berger formula or l|73l) . 



PROOF OF SEPARABILITY 

A separable expansion is always possible for local in- 
teractions provided that the potential range is finite 
(V(r) = for r > R). We follow the presentation of 
Ref. |l7j and start from the Lippmann-Schwinger equa- 

2 

tion with arbitrary energy u> (called off-shell if ui ^ ^), 
\1>) = \p)+GoVW). (76) 
We expand the wave function inside the potential range, 

i 

and multiply (|76f) by (Bj\V from the left to obtain 
2to|V|fli>& - (BjMp) +J2(B j \VGoV\B i )Z i . (78) 

i i 

Abbreviating fj,^ — (Bj\V\Bi), Q — (Bj\V\p) and intro- 
ducing the form factors 

(9j\ = (B j \V \g i )=V\B i ) (79) 

we can write (|78f) as a matrix equation 

<ft-Mft>)6 = C/ (80) 

i 

such that the expansion coefficients become 

& X njQ (81) 

i 

with tH = fiij — (gi\Go\9j)- Inserting (|8*T|) into (|77l) we 
have 



VI 



\9i)nj{9i\p) 



(82) 



which shows that the effect of the potential on the scat- 
tering wave function can be expanded in separable form. 



MODELS FOR FORM FACTORS OF 
SEPARABLE INTERACTION 

For different parameterizations of the form factor g q 
we can compute the unrestricted two-particle scattering 
expression of the T-matrix i|17fl 



which is controlled by 



XJ(P, O 



P 



2tt ! 



1 Tj 



The phase shift is then given by 



tan if 



ImT 



RcT 



AJi 



1 - A Jr. 



(83) 



(84) 



(85) 



and the parameter of the potential can be linked to the 
scattering length ao and the range of the potential ro via 
the small wave- vector expansion 



1 r 2 

p cot cp = 1 p 

ao 2 



(86) 



The T-matrix possesses a pole, the bound state energy 
Q = —Eb is given by 



1 = XJ{-E b ). 



(87) 



It will be useful to use a dimensionless scaled potential 
strength £ = mA/Air 2 f3 3 in the following. 

Exponential function 

As a first example we choose a form factor 

g p = r 2 e- p2 ^ 2 (88) 



and obtain for SI > with y' = y/mfl/2/3 2 

XJ{n) = ~iV2^ (l+F (y') + iy^y'e- y ' 
and the Fred-Conte function 

OO 2 

2z 2 f e~ x 
F(z) = — dx— -. 

V7T J X 2 - Z l 



The scattering phase shift becomes 
tan(/? = — V2irt; 



y'e-y 



T + V2<(l + F(y')) 
from which we obtain the expansion 

fl V2^ + i , V&Fg-i 2 

p cot if — —p — 1 no / . — P 



(89) 



(90) 



(91) 



(92) 
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and comparing with (|86|l we can link the potential pa- 
rameter to the scattering length and the potential range 
as 



ro = 



1 

IT 



1 



2 1 

7T 7r£ 



(93) 



For negative energies O = — Eh < we obtain with y = 

-V = V mE b/zp 2 

XJ(-E b ) = -iV2^ (l - V^ye y2 eric(y)j (94) 
such that the condition for bound states becomes 

1 + CV27T f 1 - V^n/e^erfc^)) = (95) 



from which follows that only for £ < there can be bound 
states. Rewriting it we get 



|£|V2^ = 



1 



1 — y/nye y erfc(y) 



(96) 



The right hand side is a monotonically increasing func- 
tion of y starting from 1 for y = 0. Therefore, bound 
states are only possible if 



1 



(97) 



Step function 

As a second example, the form factor is assumed to be 
a step function cut-off in real space g r = Q(R— r)/Anh 2 r 
from which we have 



1 — cospR/h 



9 P = 



For the scattering case, 57 = p 2 jm > 0, we obtain 



(98) 



V - + sin £( 0OB |-2)-i(l -cos V -f 
P P P P 



(99) 



from which the relation to the scattering length and the 
potential range follows 



1 (1 1 

2 /ll 2 



(100) 



The condition of bound states, for SI = — Eh, reads 
P / 

XJ 



^(2V2 y+ l-(e-^-2)^)(101) 



so that we can obtain bound states if 

1 _ 2y/2y + 1 - (e"^ y - 2) 2 
~ I ~ 4vV 



(102) 



and since the right hand side is a monotonically decreas- 
ing function starting from 1/3 at y — we obtain the 
condition for bound states 



£<-3. 

Yamaguchi form factor 

The traditional form factor is 

1 



(103) 



9p 



(3 2 +p 2 



For £1 = — > we obtain 
A J = £/3 2 tt 



., p 2 -f3 2 -2ip/3 



2{fi 2 +p 2 ) 2 
and the scattering phase shift becomes 
( P 2 + P 2 ) 2 p 2 -{3 2 



cot (f ■■ 



-7r£p/3 3 



(104) 



(105) 



(106) 



Expanding in p and comparing with (A3) we determine 
the potential parameters from 



a 



1 / 4 



The case of Q 



-Eh leads to 



XJ(-E b ) 



2(1 + V2y) 2 
such that we obtain bound states if 
2(1 + V2y) 2 _ 

7T 

Therefore the condition for bound states is 

7T 



(107) 



(108) 



(109) 



(110) 



(111) 



INTEGRAL OF (l50l) 

The integral 1)50(1 can be evaluated analytically. To this 
end we rewrite 1(50(1 as 



sinq/ 1 

a» — o — n ( 112 ) 

y - 2z n(y) 
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with 

n n 2 

n(y) = 1 + 2-siny + 2— (1 - cosy). (113) 

We use the representation of the principal value 1/y = 
i(l/(y — irj) + l/(y + irf)) and obtain 

l 



F= \ Q5 cl +n 2 I del Re' ' 



2izc 



-2E 



Q cos 2z+£l 



2 sin 2z 

2z 



n(2z) 



n(2z) y v -2zn'{y u )J 

-2V ^ . 

^ y u -2z n'{y v ) 

(114) 



The values y v represent the complex zeros of n(y). Set- 
ting y v — ~^2it v these zeros obey 



t&iahtv = — - 



tv 



u + n 



(115) 



Observing that l/n(2z) = A 2 Z and straightforward alge- 
bra using l|115|l leads to the final representation of the 
integral as 

^ / x ^/ ^ ^sin2z. 4 , x ^ t v 2t v +Q, 

W = n {c o,2 Z+ n—)Al + 2^— Y ^ fv . 

(116) 

The special value required for (|55|l reads 

Vt v-^ 2i„ + , . 



fi + 2t„ 



where for O > the equation (|115fl has no solution and 
therefore the last sum in (|116|) and (|117fl is absent. 
Another straightforward rewriting yields the relation 



2b 2 f dx 



2b 2 f 

Fa{ib)-F a {0) = -b+—l 



„2 j_ fr2 x 



A 2 X (118) 



which proves F^(ib) + b > Fq(0) used in the discussion 
following 1 (51)1 where b — Ry/mEb 



OTHER OCCURRING INTEGRALS 



Here we calculate the integrals appearing in chapter, 



G ni n 2 (0) 



The first integral we consider is 



(0) = dvG nin2 {v) 



^3 5g*kx(q)*k 2 (-q). 

(119) 



It is convenient here to use the wave function in momen- 
tum space. We obtain for the cavity model (0 

i \ a n sin (gift sin (x+0) 7rcos0 
Xk {q) = -4ttcR q 2R 2_ x 2 " ^S(qR-x) 



(120) 



with x = kR or for the three-dimensional wave function 



**(q) 



-cR 2 



4ttH sin (gi?) sin (x 
x(q 2 R? — x 2 ) 



-x(2ir) 3 S(qR - x) cos < 



Changing variables y — qR in (|119fl we obtain 

SRc 2 

G ni n 2 (0) — ^ 2 G*o(xi,X2) 



(121) 



(122) 



and using the exponential form factor of appendix „ we 
obtain with 7 = 2R(3 

r, 1 \ fj - v 2 /-y 2 ^ 2 sin 2/sin(xi+0i)sin(x 2 -|-(/>2) 

G o{ ^ 2 )=Jdye { y*-x 2 )(y 2 -x 2 ) 





fin 



2(x\ x 2 ) 



sin x\ sin (x<i + <p2 ) cos 0i 



+c X2 ^ 7 sin 0:2 sin (a?i + <^i) cos 02 J 
+7r 3 2 ;^(x 1 +x 2 )cos(/. 1 2 e" a: i /72 . (123) 
The remaining integral can be written as 

Xl ^ x 2 



1 



-V 2 /! 2 



sm y 



1 — 2 

x l ~ x 2 





with 

I^x 2 ) 



(y 2 -x 2 ){y 2 -x 2 2 ) 



(124) 



= dy 



2 2 

e~ y /7 siny 2 _ 1 f dtc _ t 2 1 - cos2 7 i 



2 ? 
y z — x z 



4:r; 



i — x/"f 



-ix 



(-TTlmWix/j)- f(2 1 ,x/ 1 )). (125) 



Here we used the standard integral representation of the 
complex error function for W{z) = e~ z erfc(— iz) and 
abbreviated 



dt 



e 1 cos (at) 
T~6 ' 



(126) 



This integral can be solved analytically. To this end we 
observe that 



da 



f(0,b) = -irlmW(b) (127) 
f(a,b)\ a=Q = (128) 
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and find a differential equation from i|126|) 

d 2 J{a,b) + b 2 f(a,b) = -b^c- a2 /\ 



which again is a trivial integral. 



(129) 



This linear inhomogeneous differential equation with con- 
stant coefficients can be solved straightforwardly and 
with the boundary conditions 1(1281) we obtain finally 

f(a, b) = 7re- b2 Ime- iafc erf(^ - ib). (130) 



/ dpG nin2 (p) 

The second integral used in 1)56(1 is conveniently ex- 
pressed by the wave functions in spatial representation 

dp 



(27r)3 Gni ™ 2 



(?) 



J dr 9r * kl (r/2)* k2 (-r/2)(131) 



where we use the Fourier transform of 1)121( 1 

c sin cj) ( cot x sin (kr) r < R 



*k(r) = cfce 4kr cos 



cos (kr) 



r>R ' 
(132) 



For the exponential form factor of appendix B the integrals 
in l(131|) becomes 



dp vrnpc* 

^G ni n 2 (p) = r=—(*{Xl>X2)- 



(27T) 



(133) 



The integral 1(123(1 contains a term ~ <5(xi + X2) and 
will be multiplied with ((131(1 as part of the integrand in 
(157(1 . Therefore special care is needed to evaluate 1(131(1 
for the diagonal case Xi = — X2. We split therefore 
G(x\,X2) — G'(x\ 1 X2) + G"(x\) where we use for the 
first term the radial component of ((132(1 and summarize 
in the second part the difference of the angular integral 
over the wave functions and the radial one. This differ- 
ence vanishes except for the case Xi = — X2. 

Changing variables y — r/R in ((131(1 we obtain for the 
regular part G' 

1 

G'(xi,X2) = A 2 Xi A 2 X2 J dye~ y 7 sinxiy sina^t/ 


00 

+ J dye~ v2 ~< 2 sin (x x y + cj>{) sin (x 2 y + <fe) (134) 
1 

which arc trivial integrals. The irregular part G" is now 
the difference between the angular averaging of the vec- 
torial part of the wave functions ~ k\r sinfcir and the ex- 
pression of the radial part of the wave function ~ sin 2 k\r 
since the latter is already absorbed in the regular part. 
We obtain 
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